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Abstract 

Near the horizon of a black brane in Anti-de Sitter (AdS) space 
and near the AdS boundary, the long-wavelength fluctuations of the 
metric exhibit hydrodynamic behaviour. The gauge-gravity duality 
then relates the boundary hydrodynamics for generalized gravity to 
that of gauge theories with large finite values of 't Hooft coupling. We 
discuss, for this framework, the hydrodynamics of the shear mode in 
generalized theories of gravity in d+1 dimensions. It is shown that 
the shear diffusion coefficients of the near-horizon and boundary hy- 
drodynamics are equal and can be expressed in a form that is purely 
local to the horizon. We find that the Einstein-theory relation be- 
tween the shear diffusion coefficient and the shear viscosity to entropy 
ratio is modified for generalized gravity theories: Both can be explic- 
itly written as the ratio of a pair of polarization-specific gravitational 
couplings but implicate differently polarized gravitons. Our analysis 
is restricted to the shear-mode fluctuations for simplicity and clar- 
ity; however, our methods can be applied to the hydrodynamics of all 
gravitational and matter fluctuation modes. 
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The long- wavelength fluctuations of the near-horizon metric of a black 
brane in Anti-de Sitter (AdS) space and the long- wavelength fluctua- 
tions of the metric near the AdS boundary can each be described by 
a translation-invariant (effective) thermal field theory. The equations 
of motion of either of these theories are hydrodynamic equations [I] . 
The relevant parameters — such as the temperature, shear viscos- 
ity, diffusion coefficient and entropy density — are intrinsic properties 
of the horizon and, as such, should be defined strictly in terms of the 
near-horizon metric. In spite of this apparent locality, considerable ev- 
idence has accumulated suggesting that these same horizon-specified 
quantities should be used for the boundary theory. Indeed, there are 
some concise expositions on this very point [B El O Hj; with these 
having been able to establish the boundary field theory as describing 
a viscous fluid with precisely these hydrodynamic parameters. 

The AdS boundary hydrodynamics can be related via the gauge- 
gravity duality to the hydrodynamics of strongly coupled gauge theo- 
ries [5j [3] . The latter provides an interesting theoretical framework for 
studying relativistic hydrodynamics and may explain the experimental 
results of heavy-ion collisions, non-relativistic systems exhibiting su- 
perfluidity, etc. [6j. The boundary hydrodynamics has been most ex- 
tensively studied using Einstein's theory of gravity, which corresponds 
to infinitely strong 't Hooft coupling on the gauge-theory side. As for 
applying the results to real physical systems that can be described by 
gauge theories, one really needs to know the outcomes for finite val- 
ues of 't Hooft coupling. Then, since the strong-coupling expansion 
on the gauge-theory side corresponds to the derivative expansion on 
the gravity side, such an extension actually requires knowledge of the 
results for generalized theories of gravity. More specifically, one would 
first have to calculate the higher-derivative corrections to Einstein's 
gravity in string theory and then calculate the hydrodynamics of the 
corrected theory. 

Here, we will provide a general prescription on how to perform the 
second stage for arbitrarily general gravitational corrections. First, 
we will establish a more direct connection between the near-horizon 
hydrodynamics and that of the AdS boundary. This task will be 
accomplished by showing that the thermal field-theoretic formalism 
developed for the boundary hydrodynamics can be directly translated 
into a calculation that is completely local to the horizon. Impor- 
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tantly, absolutely no information about the hydrodynamics is lost in 
this localization process. To extract the hydrodynamical parameters 
— such as the shear diffusion coefficient — at the AdS boundary, one 
solves the gravitational perturbation equations with incoming bound- 
ary conditions on the horizon and Dirichlet conditions on the outer 
boundary. However, it will be clear that one can similarly place the 
Dirichlet boundary conditions on radial shells anywhere in between 
the horizon and the boundary, including on the stretched horizon. Af- 
ter accounting for the effect of the gravitational redshift, one can see 
that the shear diffusion coefficients on all the shells are equal. 

Next, we will go on to demonstrate how the aforementioned for- 
malism can be extended to any generalized (or Einstein-corrected) 
theory of gravity. The hydrodynamic parameters of interest can be 
readily identified in terms of different components of a (generally) 
polarization-dependent gravitational coupling Kn U - Recently, we have 
proposed a method for calculating these polarization-dependent cou- 
plings that is completely local to the horizon [7]. This followed from an 
examination of the two-derivative, second-order expansion of the ac- 
tion for the gravitons (hu V ), as described in [8]. In particular, the grav- 
itational coupling can be identified on the basis that the replacement 
hu, v — > K^hfn, leads to a canonical kinetic term for the /w-polarized 
graviton. (This identification will be particularly relevant to the latter 
stages of our analysis.) It follows that any theory which is sensitive to 
the different polarizations can be expected — for instance — to alter 
the shear viscosity to entropy density ratio r//s from its "standard" 
(Einstein) value of 1/47T. This claim was put into quantitative terms 

with the following proposal [7J: ^ = j- , .2 , where the subscripts 

on the gravitational couplings denote the implicated polarizations (to 
be defined more rigorously below). For Einstein's theory or any theory 
related to Einstein's by a field redefinition, k\ v = Kg = constant (to 
be precise, one half of the d + 1-dimensional Newton's constant). 

Our particular objective will be to determine the shear diffusion 
coefficient D, for a generic gravity theory; this, by cognizance of it's 
relation to the pole structure of an appropriate thermal Green's func- 
tion. It will be shown that, in general, the product DT differs from 
the ratio rj/s; in conflict with the prediction of Einstein's theory. But 
this is not at all in conflict with the basic principles at work. As 
already stressed, each of the hydrodynamic parameters in play is sen- 
sitive to different polarizations and, therefore, will react differently to 
a non-trivial deviation from Einstein's theory. What is most signifi- 
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cant is the direct correlation between any of these parameters and the 
gravitational coupling for a very specific class of gravitons. 

We will, for illustrative purposes, be somewhat specific about the 
choice of spacetime geometry and class of perturbations; however, our 
procedure (and outcomes) can readily be repeated for many other 
interesting cases. 

Let us now proceed with the proposed analysis and consider a black 
p-brane in a d+l-dimensional (asymptotically) AdS spacetime. (Note 
that d = p + 2 > 5.) The associated metric is usually depicted in the 
Schwarzschild-like form ds 2 = -£f(r)dt 2 + + g dx 2 ) , 

where L is the AdS radius of curvature, f(r) = 1 — (r/ l /r) p+1 and 
r/j locates the black brane horizon. A simple change of coordinates, 
u = r 2 /r 2 , gives us another useful form 

r 2 T 2 rhi 2 r 2 ( P \ 

where f(u) = 1 — u~^~ . The horizon and (outer) boundary are now 
located at u = 1 and u = respectively, and T = (p + l)r/ l /4-7rL 2 can 
be identified (in any coordinate system) as the Hawking temperature. 

To study brane hydrodynamics, one expands the metric, — > 
'g^y+h^y. In accordance with the standard conventions, the coordinate 
z is singled out as the direction of propagation of the graviton on the 
brane and one of the remaining transverse directions is denoted by x 
— any of which are interchangeable by virtue of the spatial isotropy of 
the brane. Obviously, any of the x coordinates could have been picked 
instead of z for the same reason. Given these conventions, ~ 
exp[— iClt+iQz] (otherwise depending only on u), where (fi, 0, 0, Q) 
is the p + 1-momentum of the graviton. 

Under a suitable choice of gauge (namely, the radial gauge h UQ = 
for any a) , it has been shown that the non- vanishing fluctuations sepa- 
rate into three decoupled classes; with these being commonly classified 
as the scalar, shear and sound modes [9]. The latter two are of partic- 
ular interest, as the diffusion coefficient for the shear viscosity can be 
directly extracted from the pole structure of the associated correlator. 
Our current attention will be directed towards the shear channel, as 
an analogue calculation for the sound modes will yield the same basic 
outcome (albeit with some additional information) but with signifi- 
cantly more technical clutter. We will, however, briefly discuss the 
sound channel near the end of the article. 
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To determine the thermal correlator in question, one can pro- 
ceed exactly as in [lOj [ITJ E2]) where much more elaborate discus- 
sions can be found. It is first necessary to identify a gauge-invariant 
combination of the shear-mode fluctuations H± x = (—1 / gttjhtx an d 
Hzx iXl 9zz)h zx . 



Here, u = Q/2ttT and q = Q/2ttT represent a dimensionless frequency 
and wavenumber respectively. Importantly, either of these parameters 
is vanishing (although not necessarily at the same rate) in the so-called 
hydrodynamic limit. 

Restricting to the radial gauge and expanding out the Einstein 
field equations to the linear order of Z, one then schematically obtains 
Z" + A(u, q 2 f, u)Z' A- B(uj, q 2 f, u)Z = , where a prime is a derivative 
with respect to the radial coordinate u. The coefficients A and B can 
be found in, for instance, eq. (3.14) of [12J for general p (although 
with conventions differing from ours) and eq. (4.26) of [TTJ for the 
p = 3 case. What is important, for our purposes, is not necessarily 
the explicit structure of the coefficients but that q appears uniquely 
in the combination q 2 f(u). (Recall that /(it) appears in the metric 



One is instructed to solve this differential equation subject to a spe- 
cific pair of boundary conditions. Firstly, the solution is constrained 
to be an incoming plane wave at the horizon it = 1. To the linear 
order of a perturbative expansion in u and q, this condition imposes 
a solution of the form fl~2] 



Here, C is an integration constant that will be fixed by imposing the 
appropriate normalization condition. 

Secondly, there is the so-called Dirichlet boundary condition, which 
has yet to be enforced. One can impose this condition at any point 
< u* < 1; although it has become standard procedure to choose 
u* = and, thus, single out the boundary of AdS as a preferred place. 
However, this choice is not imperative. What the condition does neces- 
sitate is that Z{u) {prior to its normalization) is vanishing as u — > u*, 
which in turn imposes that u = —iq 2 f(u*)/2 — c/ eq. ([3]). As an 
immediate consequence, we see that (in spite of first appearances) uj 
and q 2 f{u*) are of the same order in the hydrodynamic limit; that 



Z = qH tx + ujH; 



zx 



(2) 



©•) 




(3) 
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is, u> 2 « q 2 f(u*). (This last observation is very important to the 
discussion that follows.) One then further requires that Z be normal- 
ized to unity at u*. For the case in hand, this can be achieved by the 
choice 

c- 1 



(4) 



At this point, let us make it clear that the hydrodynamic or zero- 
frequency limit is supposed to be put into effect. Once this limit has 
been satisfied, it becomes a straightforward exercise to show that both 
Z and its correlator (see directly below) are radial invariants. (For a 
clear demonstration of this invariance, see |13j.) 

To learn about the two-point functions of the stress tensor, it is 
sufficient to study the correlator Gzz of the gauge-invariant variable 
Z. This correlator is directly extractable from the boundary gravita- 
tional action. To elucidate, one identifies the boundary residue of the 
canonical term in the bulk action, ZZ', as the correlator up to an in- 
consequential numerical factor. Following along these lines, one is able 
to deduce that Gzz = _ bm K r , Jty <™ Z(u)Z'(u) , where K is a 

dimensional constant that depends only on the metric length scales (r^ 
and V) and the gravitational coupling. Then, taking u and u* to zero 
at the end of the calculation (and employing the Dirichlet condition to 
simplify the differentiation), one finds that Gzz = K u^^rSj^^rm 
or, by invoking J 1 <C q 2 , G Z z = — . 

One can immediately notice the pole in correlator, and, hence, the 
associated divergence. This singularity is not undesirable; rather, it 
can be viewed as a well-motivated expectation from the quasinormal- 
mode perspective of brane hydrodynamics [IT]. Moreover, by using 
the standard hydrodynamic dispersion equation to = —iDk 2 + 0(k 4 ) 
and recalling the previous scaling relations (to = 0,/2irT and q = 
Q/2irT), we can determine the diffusion coefficient as D = 1/AttT. 
The diffusion coefficient allows us, in turn, to fix the shear viscosity 
to entropy density ratio: rj/s = DT = 1/4tt in agreement with the 
"usual" (but, as discussed earlier, not necessarily universal) outcome 
for brane hydrodynamics. 

Let us now restrict our attention to near- horizon physics. After 
all, the stretched horizon — defined here as the region 1 — u* 1 - 
is a perfectly legitimate choice. So far, we have made the claim that 
the Green's function should have the same functional form at both 
the stretched horizon and boundary (and all points in between), while 
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acting as a thermal correlator for a uniquely specified hydrodynamic 
system. Let us see how this actually plays out by repeating the calcu- 
lation for Gzz but, this time, taking the near-horizon limit; that is, 
imposing u, u* — > 1 (at the very end of the calculation). As before, 
all the unwarranted zeros in the calculation nicely cancel, leaving 

r v K iq 2 f(u) 

UzZ u,u™l [w 2 - q 2 f(u)} [u + iq 2 f(u*)/2] 

lim K- — - , (5) 



u.u*— ► ! 



ILU — ■ L - L k 



where the second line follows from the hydrodynamic limit. Compar- 
ing the boundary result and ©, one can immediately see a differ- 
ence; namely, the would-be pole structure has changed by a factor of 
/(«* - !)• 

We can address this new development as follows: The frequency 
u and wavenumber q are coordinate-dependent constructs that will 
naturally experience the effect of a radially dependent gravitational 
redshift. So, how can we quantitatively discern the relative redshift 
given the sensitive (non-linear) relation between uj and q? The answer 
is remarkably simple: It is the structure of the correlator itself that 
tells us exactly how to determine the relative redshift! The hydrody- 
namic limit is, in actuality, an expansion in the ratio q 2 f(u)/io — cf, 
eq. ([3|) — so that the limiting procedure will inevitably break down 
unless q 2 f{u*)/uj < 1. Let us also observe that f(u) is monotonically 
increasing from its horizon value of /(l) = to its boundary value of 
/(0) = 1. This means that, once imposed at the stretched horizon, 
the hydrodynamic expansion at an arbitrarily larger radius is only 
ensured to persist if q 2 /u> scales as l/f(u). 

Taking u to be fixed and unaffected by the redshift, so as not to dis- 
turb the incoming (horizon) boundary condition — which, when prop- 
erly enforced, fixes a precise form for the solution at the surface of van- 
ishing / — we then have that oj(u) = uj^ and q 2 (u) = q 2 /f(u). (The 
subscript b indicates the outer boundary value of the redshifted quan- 
tities.) Now, after applying the Dirichlet condition at the stretched 
horizon to the wavefunction (|3|), one can deduce a pole structure of 
uJb = anci , accordingly, a diffusion coefficient of D = 1/AttT. In 

this way, we have finally achieved full compliance between the horizon 
and boundary calculations! 

Moreover, it should be evident that, from the current perspective, 
there is really nothing special about either the outer boundary or the 
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black brane horizon. Which is to say, the Dirichlet-imposing surface 
could have, just as well, been placed at any point between the stretched 
horizon and outer boundary without jeopardizing the form of the pole 
structure and, hence, the value of the diffusion coefficient. 

Next, let us extend considerations from Einstein gravity to a gen- 
eral theory of gravity. We will proceed to show that, for a general 
theory, the diffusion coefficient is modified in a very precise way. This 
form will be verified by two independent calculations; one of which is 
based on extracting the (modified) pole of the previously examined 
correlator and a second which considers the diffusion coefficient as a 
proportionality constant in a conservation equation for the (dissipa- 
tive) stress-energy tensor. 

By a generalized gravity theory, we have in mind a Lagrangian 
Jz? = 1 a [R + XJ^c] that allows for black brane solutions of the 

form dl]). Here, ££q represents some correction to Einstein's La- 
grangian and A is a constant "tracking" parameter. Formally speak- 
ing, the correction need not be perturbative for our framework to 
apply. However, most (if not all) interesting cases in the literature 
can be formulated as such. For further details and explanations, the 
reader can consult [7] (also, [8]). 

For a generalized gravity theory, the effective coupling can be ex- 
pressed as -Mi = 1 =F I (xfrM e afe e cd , where {a, b, c, d} € {/x, v} 

\ K l*v> \ OIX ab / 

and e a i, is the binormal vector with regard to the specified pair of polar- 
ization directions. Any binormal is antisymmetric under the exchange 
of a and b, and normalized such that e a },e ab = =p2. A =p sign is only to 
be taken as negative when one of the directions (/i, v) is timelike (this 
convention ensures the positivity of the coupling). The superscript 
(0) signifies that the calculation is always made on solution and on 
the horizon. Note that, at the order of the two-derivative expansion, 
the generalized couplings can be treated as (polarization-dependent) 
constants. 

The basic premise now goes as follows: If Einstein's gravity is 
"non-trivially" modified (meaning that the corrected theory can not 
be related to Einstein's theory by a field redefinition), then the gravi- 
tational coupling is no longer as simple as k 2 e and can be expected to 
depend on the polarization of the gravitons being probed. For exam- 
ple, calculations of the black brane entropy are known to be sensitive 
to the r — t (or u — t) polarized gravitons. Meanwhile, for the other 
hydrodynamic parameters of interest, the polarization directions will 
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depend on the particular channel being probed. More specifically, the 
scalar channel depends on the x — y (polarized) gravitons, the shear 
channel depends on the x — t and x — z gravitons, and the sound chan- 
nel depends on a particular combination of t — t, t — z, x — x, y — y, 
z — z and x — y fluctuations. (As before, z represents the direction 
of propagation on the brane, while x and y are any other transverse 
directions.) 

We are now well equipped to discuss the diffusion coefficient for a 
generalized gravity theory. Our first method is based on looking at 
the pole structure of the correlator Gzz- One might well ask as to 
how this pole would be modified for a generalized theory. To answer 
this query, let us recall the identification of the gravitational coupling, 
h[iu — ► ^ij,uhfj,u, as discussed earlier in the article. On this basis, it is 
quite natural to modify eq. ([2]) for the gauge-invariant variable Z as 
follows: 

Z = qK tx H tx + <jjk zx H zx . (6) 

One can now discover the correct scaling properties of the hydro- 
dynamical parameters in a way that resembles dimensional analysis: 
First, redefine the wavenumber and (in principle) the frequency with 
a scaling operation, second, reformulate the field equation and then 
the solution in terms of these rescaled parameters and, third, read off 
the revised pole structure. Before proceeding with the first step, let 
us fix uj (as previously explained) to preserve the incoming bound- 
ary condition at the horizon. Then, by exploiting the freedom (at 
this level of analysis) to change the normalization of Z, we obtain 
Z = q^Htx + ojH zx = q H tx + ojH zx , with q = q^-. Since the cou- 
plings can be treated as constants, the solution (|3|) remains basically 
unchanged and needs only to be rewritten in terms of the rescaled 

parameter; that is, Z ~ f(u)~~ 1 + iq %^ . Now, applying the 

L ^tx J 

Dirichlet boundary condition just like before, we can readily extract 
the diffusion coefficient 

D = Iza— (7) 

for a generalized theory of gravity. 

For a second method, our result ([7]) can be shown to agree with 
the expectations of hydrodynamics as per the following argument: In 
hydrodynamics, one can obtain D by inspecting the x component of 
the conservation equation d^T^ x = for the dissipative stress-energy 
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tensor. (The relevant stress tensor is that of the p + 1-brane theory. 
See [3] for a topical discussion.) Given the symmetries of the problem, 
this conservation equation reduces to d t T tx = d z T zx . The stress-energy 
tensor can then be expressed in terms of derivatives of the local fluid 
velocity In linear hydrodynamics, T zx = —r]d z u x whereby r/ is the 
same shear viscosity that appears in T xy = —r]d y u x . To lowest order, 
T tx = — (p + p)u x such that p is the energy density of the fluid and 
p is its pressure. Now, equilibrium thermodynamics would imply that 
T tx = —sTu x , leading to dtu x = ^d 2 u x and then DT = rj/s. In 
general, however, this is not the case. Rather, T tx = — {p + p + 5)u x 
\ u r , with x being the coefficient of heat conductance. (The correction 
5 is a purely relativistic effect that allows for an isothermal flow of 
heat in accelerated matter in the direction opposite to the acceleration 
[H1CL5].) Combining these equations, one then obtains d t u x = ^d 2 u x . 
Now, using the Kubo formula, one would find x ~ 1/ K L i n t ne same 
way that r\ ~ I/k^. So that, appropriately scaling r\ — > tj(k zx / ke) 2 

and x ~^ x( K tx/ ke) 2 , one would make D scale as D — > D - 

exactly as in eq. ((7J). 

Let us briefly point out that analogous arguments can similarly be 
applied to the sound channel; from both the thermal correlator and 
hydrodynamic perspectives. A detailed account of the sound-mode 
calculations will be reported in a separate manuscript |16j . 

Finally, the informed reader might be concerned about an apparent 
conflict between our result and that of [17] . The authors of this paper 
studied a model with the (non-trivial) correction AJzfc = ^R^uapR tlV(Tp 
and found that DT = rj/s for d+ 1 = 5 (or p = 3). We will now show 
that this agreement is actually a numerical coincidence specific to this 
particular dimensionality but is not generally true. As shown explicitly 
in our previous work [7] (where we carefully analyzed the very same 
model), 1 a = \z [l + jzd{d — 3)1 and k 2 = k% . The same 

basic procedure can also be applied to obtain . 1 .■> = , 1 .■> \l — |4dl 

(Ktx) (HE) L L J 

and, of course, k 2 zx = k e . (It should be noted that none of these 
outcomes depends upon the particular choice of radial coordinate and, 
moreover, = k 2 u for any v.) It is now straightforward to compute 
and compare the ratios of interest. Namely, from 2 = J- |^ rt ^ and 

eq. (J7|), we respectively obtain 2 — -L [\ _ j^d(d — 3)] and DT = 
■gjj [l — • Clearly, the agreement that was found between these 

two ratios is specific to the dimensionality d+1 = 5. Note that, under 
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even more general circumstances, we could not expect such agreement 
for any dimensionality. 
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